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Abstract 

The mechanism of superconductivity and magnetism and their possible interplay have recently 
been under debate in pnictides. A likely pairing mechanism includes an important role of spin 
fluctuations and can be expressed in terms of the magnetic susceptibility x- The latter is therefore 
a key quantity in the determination of both the magnetic properties of the system in the normal 
state, and of the contribution of spin fluctuations to the pairing potential. A basic ingredient 
to obtain x is the independent-electron susceptibility xo- Using LaOi_a;Fa;FeAs as a prototype 
material, in this report we present a detailed ab-initio study of Xo(q5^), as a function of doping 
and of the internal atomic positions. The resulting static XoIq^O) is consistent with both the 
observed M-point related magnetic stripe phase in the parent compound, and with the existence 
of incommensurate magnetic structures predicted by ab-initio calculations upon doping. 
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I. INTRODUCTION 



Iron pnictides^ represent a challenge in the field of superconductivity, both from the 
experimental and theoretical point of view^*^. They are magnetic metals that upon electron 
(e) or hole (h) doping transform into high temperature superconductors. The prototype of 
this new family of superconductors is LaOFeAs^, which superconducts at 26K upon partial 
substitution 14%) of O with F. Several other superconducting pnictides have since been 
discovered, with a record critical temperature up to Tc ^ 55 The parent undoped 
compound for these superconductors are anti-ferromagnetic (AFM) semi-metals, with several 
competing magnetic structures lying within a few tens of a meV-~— . This is an ideal situation 
for the presence of magnetic instabilities, with a strong possibility of spin fluctuations playing 
an important role. In fact this is conflrmed by experiments, which show that the suppression 
of the magnetic instability is a necessary step to obtain superconductivityi^*^. 

This proximity to an AFM instability leads one to draw parallels between pnictides and 
cuprates. Despite several similarities these two families of superconductors show important 
differences. One of the most important differences is related to electronic correlations; 
cuprates are well known to be strongly correlated materials, where a treatment of on site 
Coulomb interactions is essential to get the correct normal ground state. On the other hand 
in the case of pnictides it was concluded by Anisimov et al.— that the electronic structure 
of these materials is consistent with a small value of U (within LSDA+[/ scheme) indicating 
that strong correlations may not be essential to describe them. In fact, unlike in cuprates, 
standard local/semi local functional within density functional theory (DFT) can describe 
the occurrence of magnetism in pnictides, although the agreement between the calculated 
and the experimental magnetic moment and its dependence upon the choice of functional 
and atomic positions has been extensively discussed^'^. 

There is not a general consensus about the mechanism leading to superconductivity in 
these materials. It seems clear, however, that the elect ron-phonon interaction alone (without 
involvement of spins) is too weak^^ to produce such a high Tc, at least within the standard 
Migdal-Eliashberg theory. On the other hand explicit involvement of the spin degree of 
freedom leads to an enhancement of electron-phonon couplingi^i^^; therefore the role of 
phonons has not been completely ruled out. 

However, the presence of a weak electron phonon coupling at least in the classic sense 
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lead to the suggestion of several alternative mechanisms which are essentially electronic in 
nature. Among these alternatives one of the most prominent is the mechanism suggested 
by Mazin et al^; inspired by the peculiar shape of the Fermi Surface (FS)^ they suggested 
an s± superconducting order parameter, having different sign on the h and e-like sheets 
forming the FS. In this way superconductivity can be driven by a strong repulsive inter- 
band interaction (like spin fluctuations). The role of spin fluctuations in the Cooper pairing 
was proposed several years ago by Berk and Schrieffer— ; this mechanism requires a non- 
conventional (i.e. not a simple 5- wave) superconducting order parameter. The eflFective 
Hamiltonian of a system close to a magnetic instability "pairs" states at k and k + q by 
an effective interaction matrix element, ultimately related to the q- vector and frequency 
dependent susceptibility x{^^^)- Both the q- vector and frequency dependence of x 
important in determining the symmetry and possibly the anisotropy of the superconducting 
gap on the FS. Experiments, on the other hand, cannot make a deflnite statement on the 
symmetry of the order parameter^'^ making it difficult to validate/invalidate this scenario. 

It is clear that the above mentioned electronic mechanisms rely crucially on the detailed 
knowledge of the susceptibility, xio.^^)- Several calculations of the static x have been re- 
ported on pnictides. The first ab-initio calculation of the x{^^ 0) did not include matrix 
elements among Bloch states^. Later calculations2i"i2i of Xo were mostly based on the tight- 
binding fits of the minimal bands manifold around E^^ without inclusion of higher energy 
inter-band transitions. This approach suffers from the fact that, even though the main struc- 
ture of the static xi^i 0) can be determined by low energy transitions (without including 
many empty bands) no definite conclusions can be reached about the structure of xio.^^)- 
Hence an ab-initio investigation of the full dynamical susceptibility, including a proper ac- 
count of the matrix elements involving Bloch states is highly desirable and is still lacking. To 
fill this gap, in this article we present an ab-initio determination of the static and dynamical 
independent-electron susceptibility of LaOFeAs, based on electronic structure calculations 
performed within DFT. The results will be a first step towards an understanding of the ma- 
terial properties and of the possible contribution of spin fiuctuations to the superconducting 
pairing. 

The rest of the paper is organized as follows: Section [TTl discusses the methodology and 
computational details. Section [TTT] contains the results. In particular, first the static suscep- 
tibility Xo(q5 0) is discussed and then the results for the real and imaginary parts of Xo as a 
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function of frequency are presented. Finally the conclusions are given in Section HVl 



II. METHOD 



The independent electron susceptibility is defined as: 

XO (q,G,G',a;) = V /nfe " fn'k+, ^ /^^ ^ ^ ^iCq+G). /y. g-i(q+G')-r ^/j^ ^ ^ 

€nk - (^n'k+q + n[uj + IT]) ^ ' \ 

(1) 

where e^k and f^k are the one-electron energies and the corresponding Fermi functions, and 
G,G' are reciprocal lattice vectors. In this work we accurately compute Xo by performing 
the summations in Eq. ([T]) using a random sampling over the Brillouin zone (BZ). We used 
^3000 independent k— points per band, chosen according to a stochastic algorithm which 
accumulates them around the FS for bands crossing Ep. The final results are obtained by 
averaging over 40 runs each containing completely independent k-point set. This procedure 
nearly completely eliminates the numerical noise, and shows a good convergence both in 
terms of number of independent runs and number of k-points within a single run. We have 
included 65 bands in order to ensure convergence with respect to the number of empty 
bands. 

The energy bands and the matrix elements for LaOFeAs have been calculated, within 
the (spin-independent) local density approximation (LDA)^ to the exchange-correlation 
functional. Our calculations have been done using the Full-potential Linearized Augmented 
Plane Wave (FPLAPW) method. This choice is necessary because of the extreme sensitivity 
of the electronic structure (in particular, of the bands close to Ep) to the method^^. No 
numerical approximation has been made in the evaluation of matrix elements in Eq. ([1]). 
Given the huge debate about the dependence of the results on the position of the As atom 
in the unit-cel^»^, in the present work we show results both for the theoretically optimized 
( ^As = 0.638) and experimental^'' (zas = 0.6513) position of the As atom. 



III. RESULTS AND DISCUSSION 

We first look at the band structure and the FS for the parent and the doped LaOFeAs. 
The ground state calculations are performed with the tetragonal unit cell, containing two Fe 
atoms, at the experimental lattice constants a = 4. 03 A and c/a = 2.166. The energy bands 
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FIG. 1: Energy bands of LaOFeAs, for the theoretically optimized (left) and experimental position 
(right) of the As atom, along various symmetry directions in the Brillouin zone. The upper panel 
highlights the region around the Fermi level. Dashed line in the upper panel indicates the position 
of the Fermi level for the 14% doped LaOFeAs. 

of LaOFeAs calculated using both zas are shown in Fig. [H Four bands cross the Fermi level 
Ef^ having predominantly Fe d character moderately hybridized with As p. These band 
structures are in agreement with previous results^*^. 

The corresponding FS are given in Fig. [2l The FS calculated with optimized and exper- 
imental Zas for the undoped compound are quite similar with two e-like cylinders around 
the M point and two /i-like cylinders warped around the F point of the BZ. However, the 
two FS differ when it comes to the third h-like manifold which appears as a 3D structure 
on the use of optimized As positions and as a 2D cylinder when the experimental atomic 
positions are used. 

Since the parent compound becomes superconducting on electron doping, it is important 
to look at the change in the electronic structure as a function of the Fluorine content. 
The doping is treated by the means of rigid band model (REM) which is a reasonable 
approximation for LaOFeAs^^. As for the crystal structure, it is kept fixed to the undoped 
values even for the doped case. This choice is justified for the material under investigation; it 
was shown by Mazin et al.^^ that the lattice parameter and atomic positions are not sensitive 
to the doping. The FS for the 14% doped compound is shown in Fig. [2l Even though the 
undoped FS calculated using optimized and experimental zas are substantially different, the 
doped FS for the two cases are fairly similar. The h-like FS sheet around F point for x = 
quickly disappears on electron doping in both cases {zas optimized and experimental) with 
only 2D tubular structures around the T — Z (/i-like) and the M — A (e-like) lines surviving. 
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FIG. 2: Fermi surface of LaOi_a;Fa;FeAs. Top-panel: at theoretically optimized position of the As 
atom with (a) x = and (b) x = 0.14. Bottom-panel: at experimental position of the As atom 
with (c) X = and (d) x = 0.14. 




FIG. 3: Real part of x(q, 0) for undoped LaOFeAs (a) on the basal plane (i.e. = 0) and (b) on 
the top plane (i.e. Qz = ■^)' 

Most noticeably the nesting between h-like and e-like FS present in the undoped material 
becomes less prominent on doping due to the increase in the asymmetry between the e-like 
and the h-like cylinders. 

A. Static response function 

In order to elucidate how this change in FS upon doping manifests itself in the response 
of the material, in Fig. [3] and [Hare shown the response function Xo{^^ 0) for LaOi-a^F^^FeAs 
at X = and x = 0.14. Looking first at the undoped case, X = (Fig. [3](a)), two main 
features are clearly visible; first is the presence of a wide flat structure around the F and a 
relatively sharper peak around M point and second is the flatness of the response function 
in the c-axis (tetragonal axis). As expected, due to the 2D nature of this material similar 
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FIG. 4: (a) Real part of xCq^O) for 14% doped LaOFeAs and (b) Axo- the difference between the 
real part of x for the doped and undoped LaOFeAs. The results are plotted on the basal plane. 

features (structures around Z and R points in otherwise flat response function) are also 
visible in the top plane (Fig. [3](b)). These main features of the response function can be 
explained on the basis of the electronic structure of LaOFeAs and by comparison with the 
non-interacting electron gas22. The flat region near the zone center (F) is due to self-nesting 
of the FS cylinders. This scenario is similar to the 2D electron gas where the circular FS 
leads to a flat response with a sharp drop to zero at 2/ci7. In the present case, however, this 
drop is smooth, due to diflFerent sizes and irregularities in shape of the FS cylinders. The 
broad peak in the x M-point is due to the strong nesting between hole and electron like 
FS sheets. The presence of a perfect nesting, which is not realised in this material^*^, would 
give a diverging response at this point. 

Doping (see Fig. [4](a)) leads to asymmetric change in the size of the electron and the 
hole FS cylinders which manifests itself as a crater-like structure around F-point which 
continuously increases on going towards M / X points in the BZ. This monotonous increase 
of Xo with \q\ indicates the predominance of inter-band transitions. As this inter-band 
nesting becomes less perfect upon doping (worse matching of e and h FS), it produces a 
volcano-like feature around the M-point. This feature, of the maximum being shifted away 
from the M-point, is indicative of a possibility of an incommensurate spin conflguration in 
this material. Interestingly, this incommensuration has been recently predicted from two 
different ab-initio calculations^*^. The flatness of Xoi^.^^) in the remaining q— space is 
consistent with the fact that all the competing magnetic structures lie within a small energy 
range, with in plane q-vector. 

In order to understand the evolution of the features around the F and M points as 
a function of doping it is instructive to consider the contribution from different energy- 
wavevector regions to Xo- To study this we separate the contributions of different intra- and 
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FIG. 5: Intra- and inter-band contributions to Xoi^^^)^ top-panel for undoped and bottom-panel 
for 14% doped LaOFeAs. (a) and (d) are e — e , (b) and (e) are e — h and (c) and (f) are h — h 
contributions. 

inter-band electronic transitions by separating the bands as follows: the first set consists 
of the two bands forming h-like FS sheets plus all the occupied bands and the second set 
consists of the all bands forming electron-like FS plus all the empty bands. Clearly such a 
grouping of the bands disregards band crossings {e.g.^ those of the e bands with filled states 
around M). Fig. [5] shows static susceptibility resolved according to this criteria. Looking at 
various contributions around the F point it becomes clear that in both doped and undoped 
cases the low momentum susceptibility is dominated by the intra-band contribution (e — e 
and h — h). Doping leads to a strong reduction in the h — h contribution, which is expected 
since the eflFect of doping is to shrink the hole FS sheets around the F point. On the other 
hand around the M point and in the rest of the plane the inter-band e — h contributions 
dominate. Upon doping this inter-band contribution does not change substantially except 
in the very vicinity of the M point where a volcano like feature is formed. This can be easily 
understood as a consequence of the change of shape of the FS upon doping; the asymmetry 
between the e and h sheets becomes stronger resulting in a reduction of the strength of the 
nesting. 

In order to investigate how the response of LaOi-a^F^^FeAs is affected by the change in 
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atomic positions, we also calculate the Xo ? the experimental position of the As atom. 
In Fig. [6] we plot Xo^^^^ (i^ the basal plane) as a function of doping, and its change relative 
to the Xo calculated at the optimized position of the As atom. We see a sizeable difference 
between the two at zero doping: using experimental positions, the peak at the nesting vector 
is much more evident and broader, while the peak around F is depressed. We can relate this 
change to the differences in the electronic structure, the 3D bands crossing around the Z 
point become 2D, but also the radius of the FS cylinders is larger in the optimized structure 
making the nesting between e-like and h-like sheets less efficient. In the doped system, on 
the other hand, this band is completely full in both atomic configurations as a consequence 
of which in the doped case, the difference induced by change in atomic positions is not as 
significant. 

On comparing the static susceptibility calculated in the present work with the previous 
results of Mazin et al.^ one notices that the height of the peak around the M point is 
significantly suppressed on inclusion of the matrix elements. This is a consequence of the 
fact that the states near contribute to Rexo selectively, weighted by their overlap integral 
and not by unity. 

It is believed that the interacting susceptibility is a key quantity in determining the ef- 
fective electron-electron interactions leading to superconductivity. Our results show that 
the inclusion of the matrix elements influences its q dependence and in particular the rel- 
ative height of the peaks; Although our calculations only refer to the independent-electron 
susceptibility , these conclusions are likely to hold also for the interacting x- Therefore, it 
is important to include the matrix elements in order to use x for assessing the validity of 
models for the superconducting mechanism {e.g. the s± wave). 

B. Dynamical susceptibility 

An overview of the dynamical response of the system is given in Fig. [71 where we plot 
Imxo, the imaginary part of Xo (q?^)? as a function of frequency and momentum, along the 
F — M — F and F — X — F lines of the Brillouin zone, for x = 0. The most prominent 
structure is a broad peak located at frequencies between 1 and 2 eV, which behaves in a 
similar manner along the F — M — F and F — X — F directions. Near the M point we also see a 
sharp and intense peak located at about 3 eV. At this point it is also worth mentioning that 
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FIG. 6: Real part of xIq^O) ^he the basal plane (i.e. Qz = 0), calculated using the experi- 
mental position of the As atom. Top-panel: (a) for the undoped case and (b) for the 14% doped 
case. Bottom-panel: (c) difference between (a) and same calculated at the theoretically optimized 
position of the As atom and (d) difference between (b) and same calculated at the theoretically 
optimized position of the As atom. 

the flat La 4/ bands visible at around 3 eV in Fig. [1] do not result into any evident structure 
in the Imxo plot, indicating the small hybridization (thereby small matrix element) of these 
states with Fe and As states. 

At small values of q and at low frequency, the value of Imxo grows linearly as a function 
of the frequency and then drops rapidly to zero (resembling linear response susceptibility 
for the non-interacting electron gas). From Fig. [Zl(a), we see that this low frequency feature 
has two components; the position of the first peak grows with q, saturating at a frequency 
of ^0.4 eV, while the second one grows up to ^ leV finally merging into the main broad 
peak. This behaviour is the same along the F — M — F and F — X — F lines, except for 
the presence of an extra low frequency structure (cj < 0.3 eV) visible around the M-point. 
A detailed view of Imxo(q, ^) and the behaviour of Rexo{c[^(^) is given in Fig. [7](b) at X 
and M, both for the undoped and for the doped systems. Even though doping does not 
significantly change the general shape of both real and imaginary parts of Xo? iii the low 
frequency region, we notice a different behaviour of Imxo(q5^) at M (see inset of Fig. 7b). 
For X = the Imxo grows linearly as a function of cj, while it only starts to grow at a finite 
value of uj for X = 0.14. The reason for this is closely tied to the nesting function; An highly 
nested FS gives rise to linear behaviour for x = 0. While in the doped case, deterioration of 
FS nesting at M, only allows finite energy excitations, leading to a finite starting value for 



10 



3 

2.5| 
^ 21 

a 

1.5| 

1 

0.51 




1.5 



(b) 

3 



1.0 2 



V 1 



X 



MliliilL iiiii!i""-^!iliiF">iliilg 



M 





, / q=X 

j' 

i\ — 
> 


























^0.1.^ [j 

/ r 

- \ %4 


/ 








0,-5 


- K 






















( 






- ''''' /} 
" // 


- ( \ 

n 






1 


, 1 , If , 1 , 


I' 1 


/\ 1 





123 
(eV 



12:^ 



FIG. 7: (color online) (a) Imaginary part of Xo{^^^) (^V), for undoped LaOFeAs, as a function 
of q and ou. Results are plotted along F — X — F — M — F directions, (b) Plot of the real and 
imaginary parts of xo as a function of u for q corresponding to the X and M points of the BZ. 
Black and red refer to x = and x = 0.14 respectively. The inset of (b) shows the low frequency 
part of Xo at M for X = and x = 0.14 (black and red dashed respectively), and at one g-point 
point close to M corresponding to the edge of the volcano-like structure of Rexo (thin red line). 

the Imxolq, Interestingly, if we move slightly away from M (on the edges of the volcano 
structure discussed above), we recover a linear trend starting at = 0. 

The band decomposition of Imxo (q?^)? reported in Fig. [8] along the F — M — F line, 
shows that the inter-band (e — h) contributions dominate the high frequency part, and are 
also responsible for the low frequency peak around M. The two dispersive peaks at low 
frequency discussed above originate from the intra-band transitions (e — e and h — h); in 
particular, the relatively high frequency branch comes mainly from electronic bands, while 
the low frequency one from hole bands. Fig. [8] clearly shows that the h — h and e — e 
contributions are quite asymmetric; this asymmetry is clearly a consequence of the richness 
of the electronic structure of LaOFeAs near E^. 
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FIG. 8: Imaginary part of Xo(Q5^) (eV) decomposed into the intra- and inter-band contributions. 
Results are plotted along T — M — T direction, (a) h — h, ih) e — h and (c) e — e contributions. 

IV. CONCLUSION 

In summary, we report detailed calculations of the independent electron susceptibility of 
non-magnetic LaOi-a^F^^FeAs as a function of doping and of the atomic positions within the 
unit cell. Our results are based on accurate electronic structure calculations within density 
functional theory, and include matrix elements from full potential linearized augmented 
plane wave method. We account properly for Fermi surface related features through an 
accurate sampling of the Brillouin zone. 

The static susceptibility is peaked around the zone center and at the nesting vector 
(M point), due to intra-band (e — e and h — h) and inter-band transitions respectively, 
and is consistent with the observed stripe AFM ordering. However, the peak at M-point 
is not as pronounced as reported in calculations with approximate or no matrix elements. 
Upon doping, the peak at M evolves into a volcano-like structure consistent with the in- 
commensurate magnetic spiral state predicted by first principle calculations. The intra- and 
inter-band analysis of the contributions to Xo{q^^) shows an e versus h asymmetry which 
may relate to the multigap character suggested by experiments. Our results could serve 
as a first step towards the definition of an ab-initio effective electron-electron interaction, 
necessary to obtain the pairing potential in pnictides. 
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